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ABSTRACT 

This report studies the secondary pattern of a perfectly conducting offset 
main reflector being illuminated by a point feed at an arbitrary location. The 
method of analysis is based upon the application of the Fast Fourier Transform 
(FFT) to the aperture fields obtained using geometrical optics (60) and 
geometrical theory of diffraction (GTD). Key features of the present work are 
(i) the reflector surface is completely arbitrary, (ii) the incident field 
from the feed is most general with arbitrary polarization and location, and 
(ill) the edge diffraction is calculated by either UAT (Lee and Deschamps) or 
by UTD (Kouyoumjlan and Pathak). Comparison of this technique for an offset 
parabolic reflector with the Jacobi-Bessel and Fourler-Bessel techniques shows 
good agreement. Near field, far field, and scan data of a large reflector are 
presented. 
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1 . INTRODUCTION 


Reflector antennas are widely used In communication satellite systems 
because of their relatively good radiation characteristics, low cost, and 
light weight, k central problem in the analysis of a reflector antenna is 
the secondary pattern computation. As sketched in Figure 1, for an inci- 
dent field from a feed located at Pj^, the problem is to calculate the scat- 

>s * 

tered field E from a known reflector E at a far-fleld observation point r. 

Several methods exist for calculating the high-frequency asymptotic solu- 

tion of E (r), as explained below. 

(1) Physical Optics Method (PO, Fig. la) [l]-[7]. The induced 

-*-l 

current on the reflector is approximated by 2n x H . An integra- 
tion of this current over the curved reflector Z gives the far 
field E®. 

(ii) Geometrical Theory of Diffraction (GTD, Fig. lb) [8]-[10]. At a 

♦ 

tar-field observation po?nt r, the scattered field E consists of 
two terms: the reflected field on ray Pj^O*^ , and the edge 

diffracted field on ray P^O^. 

(ill) Aperture Integration Method (AI, Fig. Ic) [ 1 ] , [ 1 1 )-[ 1 3] . The 

field on aperture plane I is first calculated by tracing a 

a 

r 

reflected ray PjO ?2 using geometrical optics theory and an edge 

diffracted ray PjO^^Pj using GTD. Next, we integrate the field 

over E via FFT to obtain the scattered far-fleld 
a 

The accuracy of the above three methods is discussed below. As a 
reference, let us represent the exact solution of by a high-frequency 


asymptotic series, namely 
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Figure 1 




Three niechods for calculating scattered field E from a reflector. 
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^®(r) ~ + k“^(r) + ...], for k > « (1.1) 

Then, we may summarize the accuracy and limitations of the three methods in 
the following table: 


TABLE 1 

ACCURACY AND LIMITATIONS OF PO, AI, AND GTD 


Methods 

Accuracy 

Llmltaclou 

PO 

recover and partial Xj 

not accurate for 

AI 

wide-angle lobes 

GTD 

recover X and A, 
0 1 

predicts Infinite 
field In main 



beam direction 
(caustics) 


In this report, we will study the mai ' reflector fai"-field pattern 


using the aperture integration (AI) method cor the following reaso 

(1) Unlike Che other two methods, AI gives Che near-fleld (aperture 
field) as well as the far field. Most of today's large reflector 
measurements are done In a near field range. Thus, only AI 
provides a convenient theoretical check for the near field 
measurements. 

(11) The accuracy of AI is comparable to the popular PO. The use of 
FFT in AI makes it numerically efficient. Furthermore, as will 
be discussed later, Che present AI formulation Is most suitable 
for extension to multiple (2 or more) reflectors. 


There exists an extensive list of published literature on AI - notably. 
Silver [Ij, Kauffman and Croswell [llj, Acosta (!?}, and Hwang, Tsao, and 




Han [13]. In comparison with these prior works, the present AI antilysls may 
be considered as an extension In one or more of the following areas: 

(1) The surface of the reflector 1? completely arbitrary. It can be a 
numerically specified surface. 

(11) The edge of the reflector Is not restricted to a circular curve. 

It can be an arbitrary curve lying on an elliptical cone or 
cylinder. 

(Ill) The divergence factor of the GO field is correctly computed. 

Hence, our analysis Is not restricted to feeds located very close 
to the focal point (In which case the divergence factor Is nearly 
unity and Is ignored b;' several researchers). 

(iv) The edge diffracted field is Included In the aperture field 

calculation. Near tne Incident and reflected shadow boundaries, 
two uniform theories [14]-|15J are used so that the aperture 
field Is crntlnucus from the lit to the shadow region. 

In short, the present AI analysis represents a generalized and Improved 
version of previous work. In particular, It Is amendable for a convenient 
extension for analyzing multiple reflectors. 

The organization of this report Is as follows: The description of the 

problem Is described In Section 2. The Incident field, feed coordinates, 
and power radiated 'rrom the feed are covered in Section 3. Sections 4 and 
S give a step-by-step procedure to compute the various contributions that 
make up the aperture ^ield, namely, the geometrical optics and edga- 
dlffracted fields. Section 6 covers the aperture field theory, the applica- 
tion of the FFT, and gain normalization. Numerical results and concluding 
remarks will be presented In Sections 7 and 8. 
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2. DESCRIPTION OF PROBLEM 

The geometry of the problem under consideration is sketched in 
Figures 2 and 3. A reflector £ is illuminated by the incident field from 
a point source at P^. The problem is (a) to determine the high-frequency 
asymptotic solution of the total field at an observation point on the aper- 
ture grid Z as shovm in Figure 2, and (b) to determine the secondary 
a 

pattern using the FFT as depicted in Figure 3. In this section, we shall 
describe the various elerae..ts Involved in the problem. 

2.1. Coordinate Systems and Time Convent ion 

The main coordinate system is the rectangular system (x,y,z), whose 
origin and orientation are arbitrarily chosen. In calculating the edge 
diffracted field which involves the boundary of the refle>^tor, we employ 
a primed rectangular system (x' ,y" .z"), whose relation with (x,y,z) is 
explicitly stated later. The feed coordinate system (x^,y^,z^) required to 
describe the polarization and incident field is related to the mai ; coor- 
dinate system by Eulerian angles [16]. This will be discussed in detail in 
Section 3. The time factor is expC+jwt) and is suppressed throughout. 

2.2. Source 

We assume that the source has a well-defined "phase center" at P^ with 
coordinates (x^.y^jZ^^, and radiates a spherical wave denoted by (^^,E^). 

If the feed is an array, it is necessary to consider each element in the 
array separately and suoerimpose their final scattered fields. 

2.3. Reflector Z 


The perfectly conducting surface is described by the equation 
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Figure 2. A reflector ^ with edge F being llluainated by 
the Incident field from a point source at P^. 




Secondary 

Pattern 


1 


I 

j 



Figure 3. 


Secondary pattern using the Fast Fourier Transform. 
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2 =« f(x,y) , for a < X < b and c < y < d (2.1) 

It is not necessary to know the analytical form of the function f(x,y). In 
fact, the present computer program requires only a set of discrete data 
points (*n’^n’^n^ with n » 1,2, ...,N as the description of £. Those points 
are fitted by a cubic spline which gives automatically first and second 
partial derivatives of f, namely, 3f/3x, 3f/3y, 3^f/3x^, 3^f/3x3y, 3^f/3y^. 
There are two requirements on the cubic-spline fit: (1) the data points 

can be distributed over a random grid, but they must be dense enough to 
describe the fine details of I; (11) the domain of the data points 
(a < X < b,c < y < d) must be somewhat greater than the area defined by the 
boundary T of the reflector. Thus, we must know surface I in the shaded 
region In Figure 4 as well. Typi>_ally, the "width” of the shaded region is 
about 3 to 4 wavelengths. Our present program contains an extrapolation 
subroutine, which automatically extends Z outward "smoothly to obtain the 
necessary data points In the shaded region." The final scattered field, 
for all practical purposes. Is Independent of the surface outside F. 
Consequently, the exact manner In which the extrapolation is done Is 
unimportant. 

2.4. Boundary F 

Two types of reflector boundaries are most frequently used In prac- 
tice, and they receive our special attention. 

(a) Cylinder Case. In the first case, F Is the Intersection of sur- 
face Z and an elliptical cylinder (Figure 5a). The parameters of the 
cylinder are: 



8 
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(x ,v ) ■ center of the ellipse 
c ■ c 

(Kj,K 2 ) * semiaxis alon^ (x,y) direction. 

Using a point on the axis of the cylinder as the origin, we introduce the 

second rectangular coordinates system (x' ,y" ,z' ) such that 
■ 

X'' = X - X 

c 

\ y' = y - ( 2 . 2 ) 

z" -■» z 


The curve F may be described by a parametric equation with parameter , 
0 < ♦' < 2w, 




= p' cos 
= P' sin 

g^(^') = f(x = + gj.y = y^ + g 2 > 


(2.3) 


where 


P 



’ 

cos 

2 

' 'I 

sin 

2 




) 



(2. A) 


The projection of F on the x-y plane is always an ellipse. 

(b) Cone Case. In the second case, F is the intersection of surface 
Z and an elliptical cone (Figure 5b). The axis of the cone lies in the 
y' - z' plane and its parameters are 

(x = 0,y = 0,z = -p) = tip of cone 

0^ = inclination angle of cone axis measured from z-axis 
( 0 j, 02 ) = half-cone angles in the x' - z' and y' - 


plane 
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Using the cone tip as the origin, the primed coordinates (x^,y',z^ are 
Introduced such chat 

x" ■ X 

■ -y cos - (z - p) sin 9^ 

z' =• y sin 9^ - (z - p) cos 9^ ( 2 . 5 ) 

Using parameters the curve F is described by 
x' * ” P' cos y 

r: ■ y' * g^^') “ p' sin 

z' “ g 3 (<>') = P'[(cot 9j cos + (cot 92 sin (fr-f (2.6) 

To determine p' as a function of we must solve the following nonlinear 
equation: 

f(x,y) - p + y' sin 9^ + z' cos 9^ « 0 ( 2 . 7 ) 

For a given (()', there is a unique root p' from (2.7). The pair 

gives the desired relation p' * p'(<^''), which is fitted by spline functions. 

The projection of T on the x-y plane is. In general, a pear-shaped curve. 

(c) Arbitrary Case. In addition to the above two frequently used 
special cases, T may be an arbitrary curve described by 

x' = =■ g2(<l'').z" * gjiy) 

where (x',y',z") is related to (x,y,z) by either (2.2) or (2.5). The func- 
tions (gj.g 2 .g 3 ) can be specified either analytically or numerically. 
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2.5. Aperture Grid Points 

A general point on the aperture plane Is denoted by with coor- 
dinates (x 2 ,y 2 ,Z 2 ). Each point on E^ as shown In Figure 3 Is fed Into the 
computer for repeated calculations of the scattered field. It Is essential 
that successive observation points are adjacent, because of the following 
fact. In determining the reflection (specular) points on the reflertor, we 
make an exhaustive search only for the first observation point in a batch. 

From the second point on, we use the reflection point of the previous 
observation point as the initial guess for the current reflection point. 

It Is only when successive observation points are adjacf.nt that such an 
initial guess ensures fast convergence. In the cases that were considered, 
only one iteration was needed to obtain the reflection point for all obser- 
vation points other than the first observation point. 

2.6. Method of Solution 

-*•1 +1 

For a given incident field (H ,E ) from the source at Pp the symp- 
totic solution of the total field at point ?2 Is determined by 

Keller's geometrical theory of diffraction (GTD) [17]. Explicitly, the 
total magnetic field Is asymptotically given by 

Keller: H*^(P 2 > * + 0(k“^), k > • (2.8) 

Here the geometrical optics field, is of order k^ relative to the Inci- 
dent field and Is the dominant term. The second term Is the edge- ' 

diffracted field and is of order k It Is well-known that becomes 

Infinite and (2.8) falls if observation point ?2 Is close to the incident 
or reflected shadow boundary. In the latter case, we will use the uniform 
asymptotic theory (UAT) [14], [18]-[21], which amounts to replacing in ’ 

V 


W 




0 
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(2.8) by a new term (with capital D) : 

UAT: ^*^(P2) = H® + + 0(k“^) , k > » (2.9) 

Once is found, we calculate I*" from it by using the fact that H®, 

and are all the so-called "ray fields" which are locally plane waves. 

Once the tangential fields (f^ or H*") on Z are obtained, the Fourier 

a 

transform of these fields will essentially produce the secondary pattern. 
This will be discussed in detail in Section 6. 
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3. INCIDENT FIELD 

The definition of the incident field is given. The incident power 
from an arbitrarily polarized feed is derived so that the secondary pattern 
has the correct gain level. Finally, the computation of the incident field 
at a point on the reflector in the reflector system is described. 

3.1. Definition of Incident Field 

The surface current at the radiating aperture of the feed element may 
be expressed as 

^(x,y) - I(y ae-^'^+ yb) (3,1) 

where (a,b,<{/) are real and 

a^ + b^ - 1 . (3.2) 

By choosing (a,b,<|;), one may obtain any feed polarization. Table 2 shows 
the values of a, b, ij; corresponding to linear and circular polarizations. 


TABLE 2 

VARIOUS FEED POLARIZATIONS 

a b ij; 


linear x 

1 

0 

0 

linear y 

0 

1 

0 

RHCP 

1//2 

1//T 

90“ 

LHCP 

1//2 

HJl 

-90" 


1 

I 

I 

! 

i 


f 


! 

C 


! 


i 




I 

j 


j 

i 







15 


The radiated electric field due to 5 given by (3.1) Is 


I - 1(0 ,♦) 


(3.3) 


where 1(6,^) Is the active element pattern. The function 1(9,^) may be 
approximately expressed by 

1(0,1^) - 0Ug(9 )(ae^'*’ cos <|i + b sin iji) + (^11^(0 )(b cos ^ - ae^'^ sin (3.4) 


where 


Ug(9) ■ E-plane active pattern 
Ujj(0) ■ H-plane active pattern 

Typically, these active patterns may be approximated by (cos 9)'^, l.e., 


Ug(9) * (cos 0) 


‘lE 


U^(9) - (cos 0) 




(3.5a) 

(3.5b) 


3.2. Incident Power Radiated 

From (3.3) and (3.4), the radiated electric field Is given by 


E ~ — [Ug(0)(ae^'^ cos (> + b sin ^ )0 + Ujj(0)(b cos ^ - ae'^'*^ sin ♦)♦] 


-jkr 




J’l' 


The power radiated, assuming forward radiation only. Is 


2tt it/2 > +, 
i-0 0-0 Zo 


r2 sin 0 


where 


Z - 120ir ohms. 
0 


(3.6) 


(3.7) 


(3.8) 


Using (3.2), (3.6), and (3.8), the Incident power radiated Is 

1 * * 

^Inc ^H^h] 


(3.9) 
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For (cos 9) 


P 


Inc 


type patterns, 

60(2qg + l)(2qj^ + 1) 


(3.10) 


3.3. Determining the Incident Field on the Reflector 

The geometry of this problem Is illustrated In Figure 6. The feed 
coordinate system (x^.y^.z^) and the main r<if lector system (x,y,z) are 
related by Eulerlan angles In the most general case, these 

Cartesian systems can be aligned by threu rotations. The angles of these 
rotations are known as Eulerlan angles. Figure 7 Illustrates the Eulerlan 
angles Yj,Y 2 »T 2 * ”1^® definitions of these angles [22] are as follows. 

Angle Yj describes a counterclockwise (ccw) rotation about the z axis which 
brings the x axis to the x" axis aligned with the line of nodes (line of 
Intersection between xy and x^y^ planes), angle Y 2 defines a rotation about 
the line of nodes In a ccw sense as Indicated so that this brings the z 
axis to z^, and angle Y^ Is another rotation about the z^ axis and aligns 

A A 

the x" axis with the x^ axis In a ccw sense. Typically, x^ and z^ are 
given. Let these unit vectors be expressed by 

A AAA 

x^. « XjX + X^y + X^z (3.11) 

A AAA 

z - z X + z,y + z z (3.12) 

r I ^ 3 

Then, the Eulerlan angles are given by 






>» 


(3.15) 


ORlQlNftl. 7A’:r ^ " 
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can 


-1 


I A ~ ucait I ■ “ ' 

3 (.,X2 - Xj.2 


From [16], It can be shown that 


* ^ A m AAAm 

[x, y- Z Y - A(x y z] 


(3,16) 


where 


COSY jCOSYj-slnYjCOSY 2 ^ 10^2 slnY ^cosY j+cosY ^co 8 Y 2 SioY ^ slny^slny^ 
-cosYjSinY2”8lrtY^co8Y2C08Y2 co<jy jC08Y2COSY2-sinY jSlnY 2 3irty2CosY2 


slnY jSlnY 2 


-COSY jSlnY 2 


COSY . 


(3.17) 

Thus, a point with coordinates (x,y,z) In the main coordinate system on 

F F F 

£ has coordinates (x ,y ,z ) in the feed coordinate system given by 


F 


' 


x 


K - Xj 


F 

y 

m 

- A 

y - yi 

(3.18) 

F 




z 


' ■ 



The corresponding spherical coordinates are 

r / F.2 ^ , F,2 ^ , F,2i 1/2 
- [(x ) + (y ) + (z ) ] 


0^ - cos ^(z^/r^) 


, , -1, F, F. 

- tan (y /x ) 


(3.19) 
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From the feed function given by (3.6), the incident E-_‘ield may be 
obtained. To find the incident K-field, the following equation is used: 


♦ 

r X E 


- ["r "8 H* 1 --TT 

f f f o 


inc 


(3.20) 


where is given by (3.8). Next, the incident field is converted from 


spherical components to rectangular components using 


yf 


L Zf J 


cosOi 


-s ind , 


: 



f 

F 



f 



- 

L^fJ 


(3.21) 


Finally, the incident H-field in the ^*f*yf»2f) system is converted to the 
(x,y,z) system using the following equation 




H 

X 

*T 

^f 

H 

- A 

H 

y 


yf 

H 


H 

L z J 


L Zf J 


(3.22) 


where % is given by (3.17). 
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4. GEOMETRICAL OPTICS FIELD 

The geometrical optics field in (2.8) consists of two components: 
the incident field and the reflected field which is calculated in 
this section. It should be noted that the incident field at an observation 
point on the aperture grid is taken to be zero because the incident field 
does not contribute to the secondary pattern. 

4.1. Reflection Point 

For a given source point P^ and an observation point P 2 (Figure 2), a 
reflection point 0*^ may exist on the reflector E, and we denote its coor- 
dinates by (x,y,z»f (x,y)). The vectors ‘ 

I 

• A A A 

dj » x(x - x^) + y(y - y^) + z(f(x,y) - z^) (4.1a) 

are the connecting vectors between P^ and 0^, and 0*^ and P^, respectively. 

The condition on the reflection point is that the distance (d. d~) must 

be stationary, i.e., 

1 

i 

(d, + dj) . 0 . ^ (dj * dp - 0 (4.2) ; I 

i 

I 

i 

which is explicitly given by j 

} 

I 

{(x - x^ + [f(x,y) - - *2^ ■ ^2 1 ^ I 

^ I 

J I 

I 

I 

I 

-J- {(y - yi> + [f(x.y) ~ I iff + l<y - + [f(x,y) - z^ J = o ! 

1 2 

(A. 3a) I 

- - - - 0 
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A root (x,y,z“f) of the two nonlinear equations in (4.3a) gives the loca- 
tion of a reflection point. For a given Pj and P 2 , there may be none, one, 
or more than one reflection point. It may be shown that (4.3a) is equiva- 
lent to the satisfaction of Snell's law. 

The system of equations (4.3a) can be also satisfied if P^ 0*^, and P^ 
are collinear. Such a spurious root may be eliminated by an additional 
condition 


fx - X 


1 . jy-yi 


^2 


7 


2 - z, 


z - 2-1 


J i 


> « 


*2 J 


(4.3b) 

where 4 is a small positive number. We set 6 » 0.0001. 

A root of (4.3) may or may not fall inside the boundary T of the 
reflector (Figure 4). Thus, for each root (x,y,z=f) or its corresponding 
coordinates (x',y'’, 2 ") in the primed system, the following test must be 
performed. If 


/(x')^ + (y')^ < Aj 


(4.4) 


then the root is inside f and it is indeed a reflection point on the 
reflector. If (4.4) is not satisfied, then the root should be discarded. 
The parameter in (4.4) is given by 


‘1 O’*)'] |(^7j " , 

if r lies on an elliptical cylinder (Figure 5a); and 

1/2 


(4.5a) 


= |z'|[(x')^ + vy')^] [(x' cot 0 j)" + (y' cot 82 )“] 




- 1/2 


(4.5b) 



if r lies on an elliptical cone (Figure 5b). 
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4.2. Formula for Reflected Field 

The reflected magaetlc field at is given by 

H*^(P 2 ) - (DF) e - 2 [h^(o’^)» n] N} (4.6) 

♦ 1 

which Is given In terras of the Incident field H at the reflection po .nt 
0*^, the surface normal N of the reflector at 0*^, and a divergence fac’^or 
DF. We choose N pointing toward the source; thus, (N»z) is always gre-»tar 
than zero (Figure 8). Explicitly, N is given by 


N = A(-f X - f y + z) 

X y 


(* 4 . 7 ) 


2 2 - 1/2 

where A = +(f^ + f^ + 1) and the subscript x of f^, for example, means 
partial derivative with respect to x. The divergence factor in (4.6) is 


DF = 


/l + (d2/Rj) /l + (d^/R^) 


(4.8) 


where the square roots take positive real, negative imaginary, or zero 
value (so that DF is positive real, positive imaginary or infinite). 
(Rj.R^) are the principal radii of curvature of the reflected wavefront 
passing through 0^. Their computation is given next. 

4.3. Curvatures of Reflected W .vef ront 

We use the formulas given in Section IV of [18] for calculating 
(Rj,R 2 )» The three orthonormal base vectors of the Incident pencil are 
chosen to be (Figure 8) 

y * x^ x(z - Zj) - z(x - Xj) 

V X vi| " Tf ^2 ^ 
y * *3 I [(z - Zj ) + (x - Xj ) ] 



(4.9a) 
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X Xj 


(4.9b) 


x(x - Xj) + y(y - y^ + z(z - z^ 
[(x - Xj)2 + (y - yp2 + (z - 


(4.9c) 


where (x,y,z) are the coordinates of the reflection point 0^ . Those of the 
reflected pencil are chosen to be 


“r “1 « 

Xj * Xj^ - 2(Xj»N)N 


(4.10a) 


A J A»AA 

*2 * *2 ” 2(x2»N)N 


(4.10b) 


*3 ° T, 72~~r x2, 1/2 

l(x2 - x) + (y2 - y) + (z2 - z) ] 


(4.10c) 


^ r r ^ r 

Note that (4.10) chosen above is a left-handed system, i.e., x^ x X 2 = -x^. 
This choice, of course, does not affect the final solutions of (Rj,R 2 ). 

The three orthonormal base vectors of reflector Z at 0*^ are chosen to be 


.j, X + zf 


(4.11a) 


y + zf 
(1 + f^)^^^ 

y 


(4.11b) 


(4.11c) 


From (4.9) and (4.11) the elements 


A' > 





0 
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• m.n- 1.2,3 


(4.12) 


can be calculated with the results 


(Z - Zj) - f^(x - Xj) 




(4.13.) 


-f (x - X, ) 

y 1/ 


12 ,, r2,l/2 ,, n2 ^ / ^2,1/2 

[1 + fy] [(x - x^) + (z - Zj) ) 


(4.iJb) 


-(x-x^)(y-y^)-f^(y-y^)(z-Zj) 
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f ,^c2,1/2,, . 2 ^,, T2~ v2,2^. ^2; .2,1/2 

H+f^] {(x-Xj) (y-y^) +((z-Zj) +(x-Xj) 1 +(y-yj) (z-z^) } 


(4.13c) 


(z-z^)^-f(x-x^)^~fy(y~y^)(z-z^) 


22 


{l+fV^^{(x-x^)2(y-y^)^+l(z-zp^+(x-x^)^]2+(y-yj)^(z-Zj)2}^^^ 


(4.13d) 


A[f (x - X, ) + f (y - y, ) - (z - z,)] 


P33 "Tj" “‘‘x'* " ^1' " '•y"-' n' ‘r 


(4.13e) 


The first four elements in (4.13a) through (4.13d) form the 2x2 matrix 


Because of the particular choice in (4.10), we have ■■ P^. The 


curvature matrix of the incident pencil is 


,-l 


1 

0 


0 

.-1 


0 


(4.14) 
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The curvature matrix of reflector E at 0^ is 


A (eG - fF) A“(fE - eF) 


A^(fG - gF) A^(gE - fF)J 


(4.15) 


where 


E-l+f^, F»ff , G»l+f^ 
X X y y 


e = -Af , f = -Af , g = Af 

XX xy ® yy 


A is defined just below (4.7) 


**r 

The desired curvature matrix Q may be calculated from the following matrix 
equation 


=r 1 T -1 =7 -1 

Q*" - + 2p33[(P^)M Q (P ) 


Let us denote the four elements of Q by 


(4.16) 


Q*" = 




j ^21 *^22 


(4.17) 


Then the desired radii of curvature of the reflected wavefront at 0 are 
given by 


„r ’ r “ 2 ^ 1 ^22^ " '^^'^11 ^22^ ‘ ‘^^^n‘^22 " ^12^21^^ 


r; R2 


(4.18) 


I 

f 




Both and are real. Their signs have the following meaning: If R^ 

is positive (negative), the corresponding normal section of the reflected 
wavefront is divergent (convergent). The same convention applies to R^. 

Two final remarks about the calculation of the reflected field: 

(i) For a given and P2, there may be more than one reflection point on 
E. Then the total reflected field is the superposition of the contribu- 
tions from each reflection point. If there is no reflection point on E, 
the reflected field is zero; (ii) If the reflection point is close to the 
boundary T, we still calculate its reflected field in the usual manner. We 
shall adjust the diffracted field later by using UAT so that the total 
field if*" in this case is correct. 
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5. EDGE-DIFFRACTED FIELD 


5.1. Diffraction Points 

To calculate the edge-diffracted field, the diffraction points on the 
boundary F of the reflector must be located first. Consider a source point 
Pj at (x',y',zp and an observation point ?2 at (x£,y£,zp, with their 
coordinates given In the primed system (Figure 9). A diffraction point 
0*^ with coordinates (x',y',z") can be determined from the law of diffrac- 
tion 

“.U 

3 4 

Here t Is the unit tangent of F at o“, and 

- x'(x' - xp + y'(y- - yp +"z'(z' - zp (5.2a) 

- x"(x 2 ' - X') + y"(y' - y') + z'U' - z') (5.2b) 

From Fermat's principle, (5.1) Is equivalent to 

W ‘‘*3 '*4> ■ 0 

From (5.2), (5.3) may be written as a nonlinear equation for unknown <|)': 


f«l - . «1 " 


®2 ■ ®2 ■ ^2 

9g2 

dj d4 

V 

3^' 

[ ‘^3 ^^4 J 

3(|)' 


®3 ■ ^1 


«3 - "2 


883 

IV 


(5.4) 


where (g^,g 2 »g 3 ) are defined in (2.3) If F lies on a cylinder, and in (2.6) 


if F lies on a cone. A root of of (5.4) determines a diffraction point 
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t 


on r. Depending on the geometry, there exist examples where as many as 
four diffraction points have been found. 


5.2. Formula for Diffracted Field 

Corresponding to each diffraction point 0^, there is a contribution 
to the diffracted field in (2.8). Following Equation (5.21) of [18], the 
formula for such a contribution reads 




(5.5) 


Here g is a cylindrical wave factor 
1 


g(x) 


“■ 

f ,11 

exp 

-J X + 

^ J 


2/ 2x X ] 

The other factors used in (5.5) are explained below. 


(5.6) 


5.3. Divergence Factor 

The square root in (5.5) as usual takes positive real, negative imagi- 
nary, or zero value. is a radius of curvature of the diffracted 
wavefront passing through 0^. It may be calculated from Equation (5.11) of 
[18] which reads 


i_+ _J^ 

'^3 sin^ 3 


i- 5 - i_ 

«3 ^ <“4 


d^ Ml 


Here 3 is the angle between tangent t and d, (Figure 9). 

•4 

is given by 
sin 3 


(5.7) 


The factor sin 3 


|/l - Q^l 


(5.8a) 


0 
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where 


Pd, 




38. 


P - + 


-ci 3 g "] 

(x' - xp Ij4- + (r' - y\) ^ t (X-- xp 5 ^ 


1/2 


]'«ll 

2 

/ h 

382 

2. 

’^«3 

2 

V J 

+ 

J 

+ 

3^' 



(5.8b) 


(5.8c) 


All the derivatives In (5.8) and In the remainder of this section are eval- 
uated at the diffraction point 0*^, whose coordinates are (x' ,y' .z"* ). The 
curvature k and normal n of the curve F at 0^ can be calculated from 
Equation (13.9) of [23]. The results are 

K - ^ (5.9a) 


J_ 

PR 




8283 


“ *3*^2 


8-j8j ~ 


gjg^' - g:;g';' 




(5.9b) 


g 


1 


where the prime on g" signifies '■he oartial derivative with respect to 

n 

and 


R - [(8283' - 8382') 




<8381' - 8183') 


+ (gjg'" - 8:^8' 




2 1/2 

) 1 '^ (5.10) 


Summarizing the results in (5.7) through (5.10), we obtain the final 
expression for the divergence factor of the diffracted field 


DFD » -======- 

✓1 + (d^/Rj) 


1 

/ 1 V G 


(5.11) 
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I 


where 


^4 d4 


3 P^(sln$)2 ^1 


x' - x' x' - x' ' 
^ + £ 


’4 j 




y' - y\ y' - y\ 


l' - z" 2' - 2' 

. i+ — 


Ig;(g'jg 2 '-g 28 'r> - «3(«2*3''*3*2'^^ 


[g2(g2«3" - «3«2'^ ■ «l<*3«r ■ * 1 * 3 "^^ f 

, 


5.4. Diffraction Coefficients 

s h 

The soft and hard Keller's diffraction coefficients D and D are 
defined In Equation (5.22) of [18], namely, 


o8,h . ;; ^r 


-sec y (♦ - ± sec y ((> + 


(5.12) 


The angles and ^ are shown In Figure 10. Because of the fact that 

A A A A 

N»z > 0, It can be shown that vector t x N Is tangent to Z at 0*^, and 
points away from (not toward) £. We calculate and from the relations 


(ProJ d ). (t X N) 

cos * (-1) 

|Proj 3^1 


(5.13) 


A A 


cos ♦ “ 


(Proj d^). (t X N) 
|Proj d, I 


(5.14) 


Here Proj d^ is the projection of d 3 on the plane perpendicular to t, 
namely. 
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^ -F ^ ~ - 

Proj - (d 2 »t)t 


x'S^^ + y'S^^ + z'S33 


(5.15) 


where 


3gi 

> N * Jl. 


* (x “ x^) “■ p 


■ ~ p w 


d3p 9g2 

S32 = (y' - y'i'> - -p~ 3 ^ 


<^2^ 3g3 


S 33 = iz' - z\) - P ^ 


Similarly Proj 3^ is given by 


Proj d^ = x'S^j + y'S ^2 + z'S 


43 


(5.16) 


where 


'41 


V 3gi 

(X 2 - X ) ^ 


d4p 3g2 

■‘2 ■ ‘'I ■ ''■> 


_i- 

. I 
*■4 
■Jt 




-t 


^43 ~ ^^2 ~ ^ ^ ~ 


3g3 

~P“ 


.d , 


The normal N of the reflector at 0 is given in (4.7). For the p’*esent 
application, it is convenient to change its base vectors to those in the 
primed coordinate system. The result is 


N = x'N^ + y'N^ + z'N^ 


(5.17) 


' f 

I 


I 



where 
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N, - , N, » -Af„ . N, « A 


(5.18) 


ir boundary r lies on an elliptical cylinder; and 


^ - -Af, 


(5.19a) 


N 2 * -A(-fy cos + sin 0 j) 


N_ =• -A(f sin 0, - cos 0,) 
3 y 3 3 


(5.19b) 

(5.19c) 


if r lies on an elliptical cone. Substitution of (5.15) through (5.19) 
into (5.13) and (5.14) gives the final formulas for calculating and 


1 2 2 2 — 1/2 
cos ♦ . ISjjM, t s- SjjKjllSj, + s '2 S 33 I 


(5.20) 


222 - 1/2 

COS ♦ . + S^2«2 ^43^31^^41 ^42 ^43^ ^^’^l) 


where 


M =■ — 
1 P 




3go 3g, 


3«i 




N 


- N. 


2 P L 3 3(|»' 1 3(>'J 






- ^ 

- 1 

-9 

K 




I 


M 3 = p 


ago 3g, 

N — - N ^ 

”1 2 d^' 


The solutions of <p and are subject to the following tests: 

(i) takes the value in the range ( 0 ,tt) if 2 . 
range (ir, 2 :r) if < 0 , where 
Tj - (-Proj 13 ). N - -S 3 jNj - S 32 N 2 - S 33 N 3 


(5.22) 
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(li) 4 Cakes the value in the range (0,ir) if ^ 0, and the range 
(a,2w) if < 0, where 

= (Proj d^).N « * ^43^3 

(iii) If Che observation point ?2 is exact ly on the incident shadow 
boundary, Chen 

I - - X = 0 (5.24) 

(iv) If the observation point ?2 is exactly on the reflected shadow 
boundary, then 

♦ + - IT = 0 (5.25) 

When the observation point is exactly on the incident shadow boundary, the 
first factor in (5.12) becomes infinite, and causes computational dif- 
ficulty. A simple remedy is to shift ?2 slightly whenever (5.24) is 
satisfied. It should be remarked Chat when ?2 is near but not exactly on 
the incident shadow boundary, is large but finite. It does not cause 
any computational difficulty at the moment. Later on, the diffracted field 
in this case will be modified to become by using Che uniform asymp- • 

♦ t F 

totic theory so that the total field H is correct. Similarly, x in 
(5.12) becomes infinite if Che observation point ?2 falls exactly on the 
reflected shadow boundary. Hence, we shift P 2 slightly when (5.25) is 
satisfied. 

5.5. Spherical Components of Incident Fi el«^ 
ti- * i 

Fields H and H in (5.5) are Che two spherical components in the 
P Cl 

^ i. ^ i. ^i. d 

directions of 0^ and a of the incident magnetic field evaluated at 0 . 
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+> 




The base vectors are shown in Figure I I and may be calculated from 


*i ♦ 

6 - (a X d3>/d3 

=» (t X 3j)/(dj sin B) 


(5.26a) 

(5.26b) 


Then it may be shown that 




-1 


P sin 6 


3gj 3g2 3g3 


3 go 


‘1 “2 3^' 3 3^' 


i5.27a) 


r 

_ 


1 

3 go 


r|«i 


■ "3^ ' y? 

i 

u 

-* 


+ 


+ H, 


3go 3g, 


3g, 3go 

w - ’‘■’I 


(5.27b) 


Here (H^.H^.H^) are the components of H in the directions of (x' .y' ,z' ). 
They are given by 


Hj = H^(0‘^) , = hJ(0*^) , «3 = H^(0‘^) 


(5.28a) 


if r lies on elliptical cylinder; and 


Hj = H^(0*^) 


»frti 


H- = -cos 8. H^(0‘^) - sin 0- K^(O^) 

2 3 y 3 2 

H- = sin 6, H^(0*^) - cos 9, H^(o‘^) 

3 3 y 3 Z 


(5.28b) 


if r lies on elliptical cone. 




ORIGINAL PAGE ^ 
OF POOR QUAUi'Y 


38 






5.6. Rectangular Components of Diffracted Field 

The diffracted field calculated from (5,5) must be eventually 

AAA 

expressed In terms of rectangular components In the directions of (x,y,z), 
so that It can be conveniently superimposed with the geometrical optics 
field in Section 4, In (5.5), the two spherical base vectors (8,a) 
(Figure 11) may be calculated from 


S = (a X d^)/d^ (5.29a) 

a = (t X d^)/(d^ sin 8) (5.29b) 


The diffracted field expressed in terras of the priraed base vectors are 


^'^(P^) = x'hj + y'h^ + z'h^ 


(5.30) 


where 


h = [D^ 6 + D® a g(kd. ) (DFD) — ^ , 

n n0 na 4 sin 8 

-y2> -ar -^2> 


n = 1,2,3 


1 Pd, sin 8 
4 


2 Pd, sin i 
4 


3 g, 3 go 

(z - zp - (x - xp 


3 Pd, sin 8 
4 


^g, 3g, 

(x' ~ X') - 3^ (y' - yp 


“t ■ d" 
4 

^2 = r 

4 


yp - apz' - zp] 


^P ~ o^(x' - x')! 


'3 ■ d7 '“ 2 '"' 

4 


^p - a^(y' - yp\ 


I 

{ 

j 


\ 


I 


\ 
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The factor g is defined in (5.6), DFD in (5.11), and sin 3 in (5.8). The 
final expressions for the diffracted field in the umprimed coordinate 
system are as follows: 

h'*(P,) =* xH*^ + yH** + (5.31a) 

2 X y z 

where 

» h, , - h_ , h'* = h, (5.31b) 

X 1 y 2 z 3 

if r lies on an elliptical cylinder; and 

= h, 

X 1 

1 = -h_ cos 9- - h_ sin 0_ 

y 2 3 3 3 

= 62 sin 0^ “ h^ cos 9^ (5.31c) 

if r lies on an elliptical cone. The corresponding diffracted electric 
field E at ?2 is calculated from 

E^(P2) = 120ir d~^[H‘^(P2) x d^l (5.32) 

5.7. Detour Parameter 

The diffracted field H calculated from (5.5) is not valid when obser- 
vation point P2 is near the reflected shadow boundary which is defined by 
(5.25). To detect if ?2 is indeed so, we may calculate the so-called 
"detour parameter" of the reflected field (Section VI of [18]). 

5 = e[k(d3 + d^ - d^ - d2)]^''^ (5.33) 

Here g Is the shadow indicator of the reflected field defined by 
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e 


+1, if Is In the shadow side of the reflected field 


-I, if Is in the lit side of the reflected field 


It may be shown that 


e 



cos y ( i|> + (j)^) 


(5.34) 


The square root in (3.33) takes positive real, negative imaginary, or 
zero value. When the caustic of the reflected field falls on the reflected 
shadow boundary, C is imaginary; otherwise, ^ is always real. Detour 
parameter C in (5.33) becomes zero when ?2 is exactly on the reflected 
shadow boundary, because the diffraction point 0*^ and reflection point 
0*^ coincide (Figure 2). Following the numerical study in [24], we take 


Id - 2 (5.35) 

as on the on-set point. Thus, if )^| > 2, ?2 is considered to be away from 
the reflected shadow boundary, and the diffracted field in (5.5) is 
valid. If i^l < 2, is considered to be near the reflected shadow boun- 
dary, and we must replace by fP as stated in (2.9). 

When 5 is small, (d^ + d^) in (5.33) is nearly equal to (d^ + d 2 ). In 
many practical problems, the reflection point 0*^ and diffraction point 
0^ may not be determined with great precision. Thus, when 5 is small, a 
direct computation of C from (5.33) can have a numerical accuracy problem. 

To circumvent this possible problem, we have given below an alternative formula 
for C when its value is small: 
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cos y ((^ + 1 ^^) 


_ ✓d.tl + (d,/R,)] 
sin 6 /2k 


A + (d2/R[)/l + (d2/R2) 


cos Y ((> + ♦^) 


sin 0 /2k 


(DF)/d, 


(DDF) 


if IC 


(5.36) 


where (DF) of the reflected field is defined in (4.8), and (DDF) of the 
diffracted field in (5.11). The derivation of (5.36) is given in 
Appendix A of [14]. In all of the following computations, we use (5.33) if 
|5| > 2, and (5.36) if |5 | < 2. 


5.8. Uniform Asymptotic Theory 

We shall calculate by the UAT developed in [14], [19], [21]. The for- 
mula reads 


H°(P2) = h'*(P2) + 


F(5) - F(5) - y (1 - e) 


h’"(P2) 


where F is a Fresnel integral defined by 


F(z) = ^jW4 i’ ^-jt‘ 


dt 


(5.37) 


(5.38a) 


and 


F(z) 


2z/it 


= exp 


-j (z^ + j) 


(5.38b) 


>0|SV 

,-r 


The factor (1 - e)/2 in (5.37) is one if P 2 is In the lit region of the 
reflected field, and zero, if P 2 is in the shadow. As expected, reduces 
to when ?2 is away from the reflected shadow boundary. This is because 
of the fact that for |5 | > », 



! 

i 

‘ *• 

I 

[F(C) - F(C)) ~y (1 - e) (5.39) 

; When ?2 is near or on the reflected shadow boundary, In (5.37) is always 

finite, and compensates exactly for the discontinuity in so that the 
total field il*" in (2.9) is everywhere continuous. 
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6. SECONDARY PATTERN COMPUTATION 
6.1. Aperture Field Theory 

Refer to Figure 3. From field equivalence principles, solutions for 
the far field may be obtained knowing the tangential fields at 2^. The 
aperture plane Z is taken to be perpendicular to the z-axls. 

d 

Let us denote the tangential electric and magnetic fields at by 

S and H , respectively. The field may be determined by using vector 
d 3 

potentials [25]. However, It is more convenient to express the far field 
directly In terras of the aperture fields. Let us define the following vec- 
tor quantities: 


f(..v) . II i («.y) 
z ® 

a 

(6.1) 

g(u,v) » // S (x,y) e-^’^^“*"'’'^^^dxdy 

Z ® 
a 

where 

(6.2) 

u “ sin 0 cos ^ 

(6.3) 

V ■ sin 9 sin ^ 
k » 2tt /\ 

0,^ ■ spherical coordinates of far field point 

Since the aperture fields are tangent to z, let 
♦ * * 

(6.4) 

f - f X + f y 
X y 

^ A A 

(6.5) 

g - X + g^ y 

(6.6) 


I 

) 


I 

* 


I 

I 


I 




I 
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From each of the 3 equivalence principle formulations, three different 
expressions arise for and of the far field, namely, 

(1) using t and H 


a a 


^9 “ + fySin* + Z^cos9 (g^cos^ - g^slrvji) J 


(6.7a) 


~ “^T /r []cos9 (fyCos<t> - f^sli>|») - Z^(gySln<^ + g^cos.> )J 

(6.7b) 


(2) using 


jk 2 e ^*^*^0089 

E» » = (g„cos((i - g slty^) 

y ^ 


"9 


2xr 


(6.8a) 


ik 


(6.8b) 


(3) using E 


-jkr 


(6.9a) 


•u “ikr - 

„ ike - cos9 ^ 

“ ' 2 ^ (fyCos^ - f^sin.^) 


(6.9b) 


This method suits large apertures (in terms of X ) because the Fourier 
transforms of the aperture field, t and g, are highly peaked In the fre- 


quency domain. 
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This method Is exact If the fields at £ were known everywhere. 

a 

However, one must truncate £^ to finite dimensions in order to employ the 
FFT. 


6.2 Fast Fourier Transform 

To employ the FFT subroutine, the integrals (6.1) and (6.2) must be 
rearranged so that the form of the integral matches the definition given by 
the documentation of the FFT subroutine. For this particular 2-dimenslonal 
FFT subroutine, the function being considered Is assumed to be (.eriodic in 
X and y with period 1 in x and y. Hence, the aperture grid as shown in 
Figure 12 must be scaled accordingly. 

Many manipulations are required in order to use the FFT. These mani- 
pulations are carried out for f^^. Similarly, this can be applied to obtain 

f , g , and g . From (6.1) and (6.5) 

/A y 

f,,(u,v) - / / E (x,y) e^‘"^“’^''y^dxdy (6.10) 

X OA 

Fl Xj 

where is the x-component of ?^(x,y) in (6,1). Using the substitutions 

T - (x-x^)/(x2~Xj) (6.11) 

X ■ (y-yp/(y 2 “yi). ( 6 . 12 ) 

( 6 . 11 ) becomes 

11 _ _ 
f^(u,v) “ ^1 / / ^ax^’'i‘*’— ^’'2"*i^’ ^1 ■^X^y2"yp) • 



)2^ + v(y2-y, )X) 

dMY 


(6.13) 


where 


j 

i 

i 


I 


i 

t 


i 

I 

/ 

/ 

1 

1 
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(x 2 -Xj)(y 2 -yp e 


j|^-(uxj+vyj) 


(6.14) 
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The expression 'inder the Integral may be 

interpreted as scaled to F(a,b) where CKaO and O^bO (see Figure 

12). F(a,b) may be approximated by 


"2 "2 


FCa.W . I I 


j 2x (mX+nY ) 


n“Nj m=Mj 


(6.15) 


where C are the Fourier coefficients obtained using the FFT subroutine, 
mn ^ 

For a 32x32 FFT, 


= Mj » -15 


^2 * ^2 ~ 


From (6.13) and (6.15), 


1 1 ^2 ”2 


f (u.v) -K ! I I I C 

^ *T 


(m\+u(x2-xi))2C + (nX+v(y2-yi))xl 


0 0 n=N^ m*M^ 


mn 


dXdY 

(6.16) 


Interchanging the summation and integral signs and noting that 


1 sln(If) 

dx * e 


/ e 
0 


xc 

X 


(6.17) 


the expression for f^ is 


N 2 M 


f^(u.v) 


K I 

n=Nj m»M 


2 jr{rH-n) sln[^mX+u(x2-xi)) ]sin(-(nX+v(y 2 -yi)) ] 

I e ^ ^ * 

^ mn 
1 


ir(mX+u(x2-xi)} x(nX+v(y2-yi)) 



(6.18) 
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where 


K « (x2-xi)(y2-yi)e 


(u(xj+x 2 ) + v(yl+y2)) 


(6.19) 


In summary, Co evaluate Che Integral (6*10), we first calculate (C I of 

nn^ 


(6.15) by FFT and then calculate fj^(u,v) via (6.18). 


6.3 Gain Normalization 

The secondary pattern is usually decomposed into two orthogonal 
polarizations* Following Ludwig's definition 3 (26], Che following unitary 
polarization vectors are introduced. 


R « 0(ae^’^cos+ + bsin^) + 4(-ae^*^sin^ + bcos^) 


A A 


C * 0(ae”^'^sin^ - bcos^ ) + ^(ae •^'’cos^ + bsin^) 

Let the secondary pattern be expressed as 
-jkr 




( 6 . 20 ) 

( 6 . 21 ) 


6.S-P- [9 0, 0,] 


( 6 . 22 ) 


The reference-polarization and cross-polarization expressions of ^ are 


Reference-pol of E = (e»(R*)*) 


Cross-pol of E * (e»(C*)*) 


(6.23a) 

(6.23b) 


The second conjugate operation in (6.23) results from the change in direc- 
tion of the field after being reflected by the reflector. 

The directivity for the reference polarization is defined by 

4x |E«Rp/Z 
' o 

p 

inc 


(9.^) 


(6.24) 



similarly, the directivity for the cross polarization is defined by 


4ir |1.C|^/Z 

Or (9.^) p ^ (6.25) 

*^inc 

where is the incident power radiated from the feed, ftoting that 

e /r factor is common in both the secondary pattern (6.22) and the inci- 
dent radiated field (3.6) the directivity formulas are 


D^(6,*) 


4n |U (ae^'^cos^ + bsin^) + U (-ae^'^sin^ + bcos$)|^ 

E 1! 

Z 

o inc 


(6.26) 


0^.(0 ,♦) 


4tt |U (ae ^*^sin^ - bcos^ ) + U (ae ^*^ 008 ^ + bsin^)j 
8 ♦ 

Z P, 
o inc 


(6.27) 


Thus, for any feed polarization (a,b,iji), (6.26) and (6.27) give the 


reference and cross-polarization directivities 
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7 . NUMERICAL EXAMPLES 

The present aperture Integration method for calculating the secondary 
pattern of a reflector entails the following steps: 

(1) Input the reflector geometry, eq. (2.1), and the boundary 
description. 

(11) Input the feed pattern as described In eq. (3.4). 

(Ill) Calculate GO contribution to the aperture field. 

(Iv) Calculate edge-dlf fraction contribution to the aperture field. 

Use either UAT [14] or UTD [15] for aperture grid points near the 
shadow boundary. 

(v) Use the FFT, eq. (6.15), to obtain the far field. 

(vl) Decompose the far field pattern Into reference-pol and cross-pol 
components using eq. (6.23). 

(vll) Use eqs. (6.26) and (6.27) to obtain the reference and cross- 
polarlzatlon directivities. 

In this section, we shall present some numerical results to establ _sh 
the numerical accuracy of the present method. Near field, far field and 
scan data for a large reflector are presented. 

7.1 . Effects of Aperture Grid Size and Location on Secondary Pattern 

The aperture field theory used to determine the secondary pattern is exact 
if the tangential fields are known everywhere on the aperture plane Z . 

When employing the FFT, is truncated. To minimize the amount of com- 
puter time opent, the size of should be as small as possible while cap- 
turing almost all of the field. To this end, a study of varying the aper- 
ture grid size and location of Z^ was performed to determine their effects on 


the secondary pattern 
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The reflector used for this study Is an offset parabolic reflector 
(see Figure 13) with typical f/D^ and f/D values: 

O.AO, X- l.CO (7.1) 

P 

The corresponding values of 12^ 02* ^ 


Qj - 14.25* (7.2a) 

^2 - 64.01* (7.2b) 


H - 0.25D 


(7.2c) 


The reflector Is being Illuminated by a y^-polarlzed feed at focus with 
10 dB feed taper. Assuming (cos 0)** type patterns, eq. (3.5), the E- and 
H'-plane feed patterns are given by 


Ug( 9) « (cos 9) 

Ojj(0) - (cos 0)“*®^ 


(7.3a) 

(7.3b) 


Refer to Figure 14. Let the aperture plane Z be located a distance L 

d 

away from the focal plane. The focal plane Is located at z > f. Plots of 
the tangential components of the aperture field along two cuts of Z^ are 
obtained as the distance L Is varied. These two cuts are 


x-cut: y“H+0.5D, z*f+L 

jr-cut: x-0 ,z«f + L 

The diameter of this test reflector Is 501. Due to the choice of feed 
polarization, only the y-component of the electric field, E^, Is plotted. 
Figures 15(a)“(c) are plots of |Ey| for a x-cut for L ■ 0, 101, and 201, 




--- ~ nit, 




. I Ey I for a X'cut, L 


Figure iSa 







Figure 15b. | Ey | for 




— Total field 
GO field 








-FIELD (V/i 






respectively < Figures 16(a)-(c) are plots of jE^I for a y-cut for L ■ 0, 
lOA, and 20A, respectively. 

From Figures 15(a)-(c) and 16(a)-(c), note that the aperture field 
does not spread as L Increases. This Is so because the feed Is at the 
focus. To demonstrate that the aperture field does spread as L Increases, 
let's move the feed toward the reflector along by an amount d^ (see 
Figure 17). Figures 18(a)-(c) are plots of |Ey| for a x-cut for d^ ■ 2X 
at L - 0, lOA, and 20A, respectively. Figures 19(a)-(c) are for the y-cut. 
Figures 18(a)-(c) and 19(a)-(c) show a slight Increase of |Ey| for grid 
points near the edge of the aperture grid. Also note that the width of the 
region where the GO field Is nonzero Increased when L Is Increased. Thus, 
to minimize the field strength for a fixed size of L » 0 was chosen. 

Let the aperture grid be W by U. Choosing the criterion that the field at 
the edge of £ Is at least 20 dB below the maximum field value on £ , W « 

fl 3 

1.14D centered at the midpoint of the projected aperture was chosen. 

7.2 DBS Antenna 

A direct broadcast satellite (DBS) antenna was designed by Lee et al. 
[2> Pattern computation programs using the Jacobl-Bessel series tech- 
nique [4] were developed by Y. Rahmat-Samll. A parabolic reflector was 
used and It Is described by 

D * dish diameter * 108.148A (7.4a) 

f • focal length * 94.867A (7.4b) 

H - offset height - 16.865A (7.4c) 

The feed Is located at the focus and the primary pattern, eq. (34), Is 
described by 
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a - 1//T , b 

- 1/*T, i|; - 90“ (RHCP) 

(7.5a) 

Ug( 0) • (cos 

6)'*' 

(7.5b) 

Ug( 0) - (cos 

6)^-® 

(7.5c) 


Using this aperture Integration technique, the secondary pattern for ^ « 0° 
was computed using GO and GTD constructions. As shown In Figure 20, two 
patterns are superimposed with Rahmat-Samll's results. One pattern used GO 
fields only and in computing the other pattern, the edge-diffracted field 
was included in the aperture field calculation. The gain and sldelobe 
levels of these two patterns are tabulated versus Rahmat-Samll's results in 
Table 3 below. 


TABLE 3 

COMPARISON WITH RAHMAT-SAMII’S RESULTS 



Rahmat-Samil 

GO + GTD 

GO 

main beam 

48.28 

48.33 

48.32 

1st sldelobe 

28.42 

28.29 

28.42 

2nd sldelobe 

22.29 

22.18 

22.93 

3rd sldelobe 

18.05 

18.02 

18.12 

4th sldelobe 

14.95 

14.96 

13.40 

5th sldelobe 

12.39 

11.85 

11.14 

6th sldelobe 

10.31 

9.05 

8.41 


Due to the limited amount of computer working space, a x 48 FFT was 
used. Despite the coarse sampling, the results in Table 3 are in good 
agreement. 

The scan performance of this reflector has been studied by Hung [28]. 
The Fourier-Bessel series technique is used to compute the secondary pat- 
tern. Using the same feed, eq. (7.5), the feed was displaced 5.8 in the 






-Xp direction to produce a pattern with peak directivity at 8 ■ 3“ (see 
Figure 21). The H-plane scan pattern was also computed using the aperture 
Integration technique with and without the edge-diffracted field. Figure 
22 depicts the H-plane scan pattern using the aperture Integration and 
Fourler-Bessel series techniques. The patterns are In good agreement. 

7.3. TRW Antenna 

A dual reflector antenna system was designed by TRW for NASA-Lewls 
Research Center. The dual reflector Is an offset Cassegrain reflector. 

The main reflector Is parabolic and it Is described by 

D “ dish diameter =• 257. 89A (7.6a) 

f » focal length » 318. 74X (7.6b) 

H » offset length » 135. 51X (7.6c) 

The reflector Is being illuminated by a y^-polarized feed at focus with 
18 dB edge taper. 

Two cuts of Z at the focal plane were taken, 
d 

x-cutr yaH+0.5D, z*f 
y-cut: X ■ 0 , z » f. 

The magnitude and phase of the y-coraponent of the electric field, are 
plotted in Figures 2i(a)-(b) for the x-cut and in Figures 24(a)-(b) for the 
y-cut. 

The secondary pattern for * 0® of this main reflector is shown in 
Figure 25. The key features of the ref erence-pol directivity plot are 
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Figure 25. Reference-pol. directivity plot 
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‘ tv 

-t 
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t: 


GO + GTD GO 


Gain (dB) 

56.95 

56.95 

HPBU (deg) 

0.276 

0.277 

SLL (dB) 

-32.00 

-32.26 


Figure Tj is a plot o‘ Che cross-pol directivity. 


The hyperboloid subref lector has a aagniflcation factor of 2. In 
approximating the performance of this dual reflector system, we employ the 
equivalent paraboloid [29]. The equivalent paraboloid is described by 

D - 257.89 (7.7a) 

f - 2(318.74) » 637.48 (7.7b) 

H » 135.51 (7.7c) 

The corresponding secondary pattern for * 0* is shown in Figure 27. 


and side lobe 

level (SLL) 

are 


GO + GTD 

GO 

Gain (dB) 

56.87 

56.88 

HPBW (deg) 

0.279 

0.280 

SLL (dB) 

-33.87 

-35.15 


Figure 28 is a plot of the cross-pol directivity. 

Scan performances of Che TRW main and :r juivalent reflector were 
studied. The scan plane chosen Is the plane and E- and H-plane scans 

were performed (see Figure 21). 

For E-plane scan, the coordinates of the feed are 


Ms* 
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X • 0 

y » -f tan (N«HPBW) cos 0^ 

2 - f - f tan (N*HPBW) sin 9^ 

where 


(7.8a) 

(7.8b) 

(7.8c) 


f « fecal length 

„ „ . -1 + 0.5Dx 

0^-2 tan ( 2 f ' 


(7.9) 


N ■> number of beamwldth's scar. 

HPBW » beamwidth of secondary pattern ( * 90®) with feed at focus. 
For H-plane scan, the feed coordinates are given by 

X - -f tan (N.KPBW) (7.10a) 

y - 0 (7.10b) 

2 - f (7.10c) 


where 


f » focal length 
N » number of beamwidth' s scan 

HPBW » beamwidth of secondary pattern ( (fr ■ 0®) with feed at focus. 

As the feed moves away from the focus, the secondary pattern degrades. 
The degradation may be characterized by peak gain loss and half-power beam- 
width. These 2 flgures-of-merit are plotted versus number of beamwldths 
scanned for the following four cases: 

(1) TRW main reflector. E-plane scan (see Figures 29-30) 

(2) TRW equivalent reflector. E-plane scan (see Figures 31-32) 

(3) TRW main reflector, H-plane scan (see Figures 33-34) 

(4) TRW equivalent reflector, H-plane scan (see Figures 35-36). 

As expected, the scan performance of the equivalent reflector is much 
better than for the main reflector due to its larger f/D value. For ti 
same number of beamwldths scanned, the peak gain loss and amount of beam 





PEAK GAIN LOSS (dB) 



BEAMWIDTHS SCANNED 


Figure 29. Peak gain loss versus beamwldths scanned for E-plane scan of TRW’ 
main reflector. 
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Figure 35, Pea^ gain loss versus beamwidths scanned for H-plane sc;?n of TRW 
equivalent reflector. 
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Figure 36. Half 'power beaowldch versus beamwldchs scanned for H-plane scan 
of TRW equivalent reflector. 
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broadealng were much less for the equivalent reflector for both and H- 
plane scans. Also, the cross-pol directivity was much lower for the equiva- 
lent reflector. Figures 37(a)-(b) and 38(a)-(b) are the reference-pol and 
cross-pol directivity plots for 6 beamwidth scan in the H-plane for the 
TRW main and equivalent reflectors, respectively. With respect to Che peak 
directivity, the maximum cross-pol directivity value is -39.15 dB for the 
equivalent reflector and only -27.95 dB for the main reflector. 


I 




9 (degrees) 


Figure 37a, Reference-pol. directivity plot of TRW main reflector for 6 
beaoivldth scan in the H-plane, 
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Figure 37b. Cross-pol. directivity plot c 
width scan in the H-plane. 
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TRW main reflector for 6 beaor* 

















8 . CONCLUDING REMARKS 


We have developed a method of calculating the secondary pattern of an 
arbitrarily shaped reflector illuminated by a feed with arbitrary polariza- 
tion. An edge-diffracted field was added to the geometrical optics field in 
the aperture field calculation. By employing the FFT, the secondary pat- 
tern is computed very efficiently. The results for the secondary pattern 
are in good agreement with those obtained by the physical optics Integral. 
Furthermore, this method can be conveniently extended to secondary pattern 
computation of multlpl'? reflector systems, which will be done in the next 
phase of this project. 


V , 
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APPENDIX A 


INSTRUCTIONS TO RUN COMPUTER PROGRAMS 
To run the computer programs on the University of Illinois CYBER 175 
system to compute the secondary pattern of an offset parabolic reflector, the 
user should perform the following 

(1) create input file TAPE2 

(2) run program BALI 

(3) run programs BPT, BFFT 

(4) for plot, use BPL0T3 

The computer commands corresponding to steps (l)-(4) are 

(first create TAPE2) 

R 

BALI 

R 

P.LOAD(BPT,BFFT); EXECUTE 
R 

BPL0T3 

The secondary pattern is stored in TAPE13. 

The input file TAPE2 comiists of 14 lines. These inputs are described 
below. 

1. H,D,f,Zo — offset parabolic reflector parameters as depicted in 
Figure 13 

2. 1 — cylinder boundary or 2 — cone boundary 

3. x^, y^, Kp K 2 — cylinder boundary parameters or P, 0p02,0|j — 
cone boundary parameters as shown in Figure 5 

YjfY 2 »Y 3 — Eulerlan angles (see Figure 7) 
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5. — source point coordinates In main reflector coordinate 
system. 

6. a,b,i|; — polarization of feed 

7. 9g*9jj — feed taper parameter for E“ .• .u h " ane primary patterns 

8. L — location of aperture plane beyond the fo^al plane as shown 
In Figure I4a (taken to be zero) 

9. type In 1 or 2 

10. XpX 2 »N X, ,x«,y, ,y^ are bounds of E (see Figure 12a), N Is the 

, , „ number of points In sampling. 

»y2 

12. 0 — use UAT or 1 — use UTD, 0 — compute GO and edge-diffracted 
field or 1 — compute GO field only 

13. 6p02,NFFP — For a theta cut, 0^ is the first value and ©2 is the 

last value of theta. NFFP is the number of 0 values for the 

secondary pattern 

14. 1 — use E-flelds only for the aperture field calculation or 2 — 
use H-fields only, PHI ~ this is the constant ^ value when 
computing the secondary pattern. 

For an arbitrarily shaped reflector, the user should use the binary ver- 
sion of MAIN Instead of BALI. However, lines 97 to 107 must be replaced with 

either the coordinates of the points that define the reflector or the new 
equations that describes the refl* .^tor. For an arbitrary reflector, line 1 of 
input file TAPE2 has no meaning. Thus, delete line 41 of MAIN and the input 
file is lines 2 through 14 of TAPE2 described above. 
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